Dynamic light scattering measurements in the activated regime of dense colloidal hard 

spheres 
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We use dynamic light scattering and numerical simulations to study the approach to equilibrium 
and the equilibrium dynamics of systems of colloidal hard spheres over a broad range of density, 
from dilute systems up to very concentrated suspensions undergoing glassy dynamics. We discuss 
several experimental issues (sedimentation, thermal control, non-equilibrium aging effects, dynamic 
heterogeneity) arising when very large relaxation times are measured. When analyzed over more 
than seven decades in time, we find that the equilibrium relaxation time, Tq, of our system is 
described by the algebraic divergence predicted by mode-coupling theory over a window of about 
three decades. At higher density, increases exponentially with distance to a critical volume 
fraction ipo which is much larger than the mode-coupling singularity. This is reminiscent of the 
behavior of molecular glass-formers in the activated regime. We compare these results to previous 
work, carefully discussing crystallization and size polydispersity effects. Our results suggest the 
absence of a genuine algebraic divergence of in colloidal hard spheres. 

PACS numbers: 05.10.-a, 05.20.Jj, 64.70.P- 



I. INTRODUCTION 

Colloidal particles arc increasingly popular as model 
systems for investigating the behavior of atomic and 
molecular materials [l[. As the typical size of a col- 
loidal particle is comparable to the wavelength of the 
visible light, one can use relatively simple, yet powerful, 
techniques, such as optical and confocal microscopy 01 
or static and dynamic light scattering to fully char- 
acterize their structural and dynamical properties. In 
contrast to molecular systems, the interaction between 
colloids can easily be tuned from repulsive (e.g. due to 
electrostatic interactions), to hard sphere- like mod- 
erately attractive (e.g. due to depletion forces Q), or 
even strongly attractive (e.g. due to van der Waals at- 
tractions). Furthermore, colloidal particles can now be 
synthesized in a variety of non-spherical shapes, and their 
surface can have non-uniform physical or chemical prop- 
erties, opening the way to anisotropic systems with di- 
rectional interactions that can be precisely tailored [1]. 

This variety of morphologies and surface properties, 
and the possibihty to study colloidal systems with rela- 
tive ease, comes with a price. Unlike atoms, colloids have 
to be synthesized and the outcome is usually a polydis- 
perse assembly of particles with different sizes. Depend- 
ing on the phenomenon of interest, polydispersity might 
be an important control parameter, e.g. when the precise 
location of phase transitions must be determined. More- 
over, the very same feature that makes accessible their 
typical time and spatial scales — their comparatively large 
size — also makes colloidal particles prone to sedimenta- 
tion. The density of a colloid typically differs from that 
of the solvent in which it is suspended. Therefore, col- 



loids experience a buoyant force, proportional to Ap, the 
difference between the density of the particle and that of 
the solvent. Sedimentation (or creaming, if Ap < 0) is 
usually not a severe issue for isolated, submicron col- 
loids. The situation is however radically different for 
colloids that form solid structures, such as aggregates, 
crystallites, or glasses, because the gravitational stress 
due to a large number of particles may add up. These 
effects can be mitigated only partially by matching the 
solvent density to that of the particles 0, [1| , since a per- 
fect matching can not be achieved in practice. Finally, a 
large particle size implies that microscopic motion occurs 
on a timescale typically much larger than for molecular 
systems. Although clearly an advantage when single par- 
ticle motion is investigated via direct visualization tech- 
niques, this is not necessarily so when slow, collective 
relaxation phenomena are studied, which can occur on 
timescales of several days. 

In this paper, we discuss some of these issues in re- 
lation with our investigations by dynamic light scatter- 
ing (DLS) of dense suspensions of colloidal hard spheres 
approaching the colloidal glass transition For hard 
spheres at thermal equilibrium, several distinct glass 
transition scenarios have been described. In a first line 
of research, the viscosity or, equivalently, the timescale 
for structural relaxation, Tq,((^), is believed to diverge al- 
gebraically: 
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This result is both predicted [l3| by mode coupling the- 
ory (MCT), and supported by Hght scattering data [nj. 
Packing fractions (pc « 0.57 — 0.59 are the most often 
quoted values for the location of this colloidal glass tran- 
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sition. A truly non-ergodic state is also often reported 
for larger <p [111, [11, [Mill . 

Several alternative scenarios [H, [H, [13, [Il| suggest a 
stronger divergence: 
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Equation ([2]) with (5 = 1 is frequently used to account 
for viscosity data [H, [l^j because it resembles the Vogel- 
Fulcher-Tammann (VFT}form used to fit the viscosity of 
molecular glass- formers |2l| . with temperature replaced 
by ip. Moreover, it is theoretically expected on the basis 
of free volume arguments [l^, which lead to the iden- 
tification = f^rcp, the random close packing fraction 
where osmotic pressure diverges. Kinetic arrest must oc- 
cur at (picp (possibly with 5 = 2 [IBl), because all parti- 
cles block each other at that density [l^, [H, [2^ . This 
is analo gou s to a T = glass transition for molecular 
systems [24| . 

Entropy-based theories and replica calculations [13, [13 
predict instead a divergence of Tq, at an ideal glass tran- 
sition at Lpo < (/^rcp, where the configurational entropy 
vanishes but the pressure is still finite. This is anal- 
ogous to a finite temperature ideal glass transition in 
molecular glass- formers 2j|. In this context, the connec- 
tion to dynamical properties is made through nucieation 
arguments [25[ yielding Eq. with 6 not necessarily 
equal to unity [26[. Here, the relaxation time docs not 
diverge because particles are in close contact, but be- 
cause the configurational entropy counting the number 
of metastable states vanishes. 

In molecular glass-formers where dynamical slowing 
down can be followed over as many as 15 decades, the 
transition from an MCT regime, Eq. ([T|). to an activated 
one, Eq. (21); has been experimentally demonstrated [2l[. 
For colloidal hard spheres, the situation remains con- 
troversial, because dynamic data are available over a 
much smaller range 

[ill, [13, [H, [13, typically five decades 
or less. Crucially, only equilibrium measurements for 
(p < ^Pc were reported, leaving unknown the precise na- 
ture and location of the divergence. A confusing feature 
of the colloidal glass transition, therefore, is that the 'ex- 
perimental glass transition', which we denote pg, and 
which occurs when equilibrium relaxation times become 
too large to be confidently measured experimentally, oc- 
curs near the fitted MCT singularity, Lpg ^ (fc, while the 
two crossovers are well-separated in molecular systems. 

At the theoretical level, there is also an on-going debate 
about the nature of the glass transition in colloidal hard 
spheres systems, theoretical claims exist that the cutoff 
mechanism suppressing the MCT divergence in molecular 
systems is inefficient in colloids due to the Brownian na- 
ture of the microscopic dynamics, suggesting that MCT 
could be virtually exact for colloids |28j . This vi ewp oint 
is challenged by more recent MCT calculations [U, [s^l , 
and by computer studies of simple model systems where 
the MCT transition is similarly avoided both for stochas- 
tic and Newtonian dynamics [sil . [3^ . [ssj , directly empha- 



sizing that the cutoff mechanism of the MCT transition 
in molecular systems is different from the one described 
in Ref . [1^ , and is likely very similar for hard sphere col- 
loids and molecules. The only physical ingredient missing 
in these theoretical works is the inclusion of hydrody- 
namic interactions, which are always supposed to play 
little role at very large volume fractions. It would be 
very surprising if hydrodynamic interactions in colloids 
could suppress activated processes and make the MCT 
predictions exact. 

In this article we discuss in detail all these issues. 
In a short version of this work [^, we claimed that we 
had been able to detect ergodic dynamics for a colloidal 
hard sphere system at volume fractions above the mode- 
coupling transition, ipc, because our experiment was able 
to cover an unprecedented range of variation of the equi- 
librium relaxation time. In the present article, we discuss 
in detail the several challenges we had to face in order 
to obtain these data, and we argue that our experimen- 
tal results should apply quite generally to colloidal hard 
spheres. We first describe our sample and experimen- 
tal set-up in Sec. [TTl In Sec. IIIII we specifically discuss 
the difficulties associated with measurements of long re- 
laxation timescales. In Sec. llVl we analyze our results for 
the equilibrium dynamics, and compare them with previ- 
ous experiments and a new set of numerical simulations. 
Finally, we conclude the paper in Sec. El 



II. COLLOIDAL PARTICLES AND 
EXPERIMENTAL SETUP 

A. Sample preparation 

The particles are poly(mcthylmethacrylate) (PMMA) 
spheres of average radius R = 130 nm, sterically sta- 
bilized by a thin layer of poly(12-hydroxystearic) acid 
of thickness ~ 10 nm. They are physically similar to 
those used in previous studies of the glass transition (see 
e.g. [3, [ill, [3)- A major advantage is that they are 
slightly smaller than in previous studies. This implies 
that relaxation in the dilute limit is faster, and that we 
can cover a broader range of volume fractions. The par- 
ticles are polydisperse (standard deviation of the size 
distribution normalized by the average size a « 10%, 
as obtained by a cumulant analysis of DLS data [s^). 
in order to prevent crystallization for at least several 
months, much longer than the typical experiment du- 
ration. They arc suspended in a mixture of cis/trans- 
decalin and tetralin (66/34 w/w), whose refractive index 
riD = 1.503 at r = 27 °C closely matches that of the col- 
loids, thereby minimizing van der Waals attractions and 
allowing light scattering experiments to be performed in 
the single scattering regime. Prior to each measurement, 
the suspensions are vortexed for about 6 h, centrifuged 
for a few minutes to remove air bubbles and then kept on 
a gently tumbling wheel for at least 12 h. The waiting 
time or sample age, iw, is measured from the end of the 
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tumbling. When several measurements are repeated on 
a sample at a given volume fraction, the dynamics arc 
re-initialized by tumbling it on the wheel overnight, with 
no further vortexing or centrifugation. 

Samples at various volume fractions are prepared as 
follows. A stock suspension is centrifuged in a cylindri- 
cal cell for about 24 h at 1000 g, where g is the accelera- 
tion of gravity, to obtain a random close packing (RCP) 
sediment, whose volume fraction, Lp, is estimated to be 
f = </?rcp ~ 0.67 according to numerical results [1^. It is 
crucial to remark that this value of (^rcp is affected by a 
large uncertainty, since it depends on the details of the 
particle size distribution, but also because the polymer 
layer covering the particles may be slightly compressed 
during strong centrifugation, implying that volume frac- 
tion is not accurately known at this stage of the prepa- 
ration. 

Solvent is then added to, or removed from, the clear 
supernatant in order to adjust the overall volume frac- 
tion io if K 0.4. The suspension thus obtained is used as 
a mother batch from which individual samples are pre- 
pared at the desired volume fraction by adding or re- 
moving a known amount of solvent. All volume fractions 
comparative to that of the initial batch are obtained with 
a relative accuracy of 10"'*, using an analytical balance 
and literature values of the particle and solvent densities 
(1190 kgm-3 and 930 kgm'^, respectively) [l^. 

Although the relative volume fractions within our ex- 
perimental data arc known with high accuracy, the abso- 
lute scale of ip is far more difficult to estimate precisely, 
as there is no direct way to measure Lp which is not prone 
to uncertainty. We discuss this issue in more detail in 
Sec. HID] below. 



B. Dynamic light scattering setup 

We use dynamic light scattering to probe the dynam- 
ics of our concentrated colloidal hard spheres. In a DLS 
experiment one measures the autocorrelation function of 
the temporal fluctuations of the intensity scattered at an 
angle 9 Q. This allows the dynamics to be probed on a 
length scale ~ 27r/(7, where q ~ 2ko sin(0/2) is the magni- 
tude of the scattering vector, with fco the wave vector (in 
the solvent) of the incident light, usually a laser beam. 

Dynamic light scattering experiments are performed 
on samples thermostated at T = 27±0.1 °C, using 
both a commercial goniometer and hardware correlator 
(Brookhaven BI9000-AT) , to access the dynamics on time 
scales shorter than 10 s, and a home-built, CCD-based 
apparatus to measure slower dynamics. The CCD appa- 
ratus is described in Ref. [36j . and so we simply recall 
here its main features. The source is a solid state laser 
with in vacuo wavelength 532.5 nm and maximum power 
150 mW. The sample is held in a temperature-controlled 
copper cylinder, with small apertures to let the incom- 
ing beam in and out and to collect the scattered light. 
The beam is coUimated, with a diameter of 0.8 mm. All 



data reported here are obtained at a scattering angle 
6* = 90 degrees, corresponding to g = 25.08 /xm~* and 
qR = 3.25, slightly smaller than the interparticle peak in 
the static structure factor. The collection optics images 
a cylindrical portion of the scattering volume onto the 
CCD detector, the diameter and the length of the cylin- 
der being approximately 0.8 and 2 mm, respectively. 

We process the CCD data using the Time Resolved 
Correlation (TRC) method [IJ , which allows us to char- 
acterize precisely both equilibrium and time-varying dy- 
namics. The degree of correlation, c/, between pairs of 
images of the speckle pattern scattered by the sample at 
time tw and + t is calculated according to 

, _ (/p(tw)/p(tw+r))^ _ 
''^*^'"^^(/p(M>,(/p(iw + r))^ 

where Ip is the intensity at pixel p and the average is 
over all CCD pixels. Data are corrected for the un- 
even distribution of the scattered intensity on the de- 
tector as explained in Ref. [11]. For non-stationary dy- 
namics (e.g. during sample aging) the two-time inten- 
sity correlation function, g2(t^,t^ -I- t) — 1, is obtained 
by averaging ci{t^,T) over t^, choosing a time win- 
dow short enough for the dynamics not to evolve sig- 
nificantly. For age-independent dynamics (e.g. when 
equilibrium is reached), time invariance is fulfilled and 
<?2(^w,iw + t) — 1 reduces to the usual intensity correla- 
tion function g2iT) — 1. 

In order to obtain the full correlation function at all 
relevant time delays, we first use the CCD apparatus to 
measure the dynamics for time delays r > 0.1 s. By ap- 
plying the TRC method, we monitor the evolution of the 
dynamics, until an age-independent state is reached. The 
sample is then transferred to the goniometer setup, where 
the fast dynamics (r < 10 s) are measured using a point- 
like detector and averaging the correlation function over 
time. The intensity correlation functions measured in the 
two setups arc then merged by multiplying the CCD data 
by a constant, so that the two sets of data overlap in the 
range 0.1s<r<10s. Finally, the full intensity correla- 
tion function is rescaled so that 52(t) — 1 ^ 1 for r 0. 
Although in this paper we focus on the slow dynamics 
and hence show only the CCD data, we point out that 
the goniometer measurements are still necessary in order 
to properly normalize the intensity correlation function. 
For (p > 0.55; we find that the dynamics are too slow 
to obtain properly averaged data when only a time av- 
erage is performed. We thus adopt the so-called brute 
force method (39| for the goniometer measurements: a 
large number of intensity correlation functions are col- 
lected (typically 200 to 500), the sample being turned in 
between two measurements in order to probe statistically 
independent configurations. Additionally, the collection 
optics are designed in such a way that about 10 coherence 
areas, or speckles, are collected by the detector, further 
enhancing the statistics. 
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C. Coherent vs. incoherent scattering 

The intensity correlation function is related to 
f{tv!,tvj+T), the (two-time) intermediate scattering func- 
tion (ISF) or dynamic structure factor, by the Siegert re- 
lation [401 : / = P~'^y/g2 — 1; where /3 < 1 is a constant 
that depends on the collection optics. The ISF quantifies 
the evolution of the particle position via 



^ N N 



,iq-[r,(t„+r)-rfc(t„)] 



(4) 



3 = 1 fc=l 



Here, the double sum is over N particles, rj{t) is the po- 
sition of particle j at time t, q is the scattering vector and 
brackets indicate the same average as for 52 • Note that 
in general DLS measurements provide information on the 
collective (or coherent) ISF, as indicated by the double 
sum over j and k in Eq. For optically polydisperse 
samples, however, the refractive index of the solvent may 
be tuned so that the ISF reduces to its self (or incoher- 
ent) part, i.e. only terms with j = k contribute to the 
r.h.s. of Eq. (g]). Our sample is optically polydisperse, 
because the PMMA core is polydisperse in size and the 
grafted PHSA layer has a refractive index that is differ- 
ent from that of the PMMA core, so that the average 
refractive index of a particle varies with its size. Since 
the size polydispersity is only moderate, one can safely 
assume that the scattering power of a given particle is 
independent of the particle position. Under this assump- 
tions, the ISF can be written as the sum of a "full" and 
a "self" (incoherent) term [s^l- The contribution of the 
full term vanishes if one matches the solvent refractive in- 
dex to the mean refractive index of the particles. In the 
experiments performed to calibrate the absolute volume 
fraction by measuring the ip dependence of the short-time 
dynamics for (p < 0.2 (see Sec. HID I and Fig. [1] below), 
we carefully adjust the refractive index of the solvent so 
as to fulfill this condition. Additionally, these measure- 
ments are performed at qR = 4.0, for which S{q) « 1, 
where the only contribution to the ISF is the self part 
even if the average index is not matched (3^ . Thus, in 
these runs we probe uniquely the self part of the ISF, as 
required for applying the calibration method described 
below. 

In all the other measurements reported in this pa- 
per and in Ref. we deliberately increase by a small 
amount the optical mismatch between the particles and 
the solvent, in order to increase the intensity signal to be 
detected by the CCD camera. By comparing the mean 
intensity for the slightly mismatched sample to that mea- 
sured for the best match, we find that the relative weight 
of the "full" and "self" terms of the ISF are 35% and 
65% respectively. The "full" term itself is a combina- 
tion of the self and distinct parts of the ISF, the relative 
weight of the self part being of order 1/S{q) [s^l- At the 
scattering vector where the measurements are performed 
{qR = 3.25), S{q) < 2 for all studied volume fractions. 
Therefore, the relative weight of the self term in the mea- 
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FIG. 1: (Color online) Volume fraction dependence of the 
dimensionless short time self-diffusion coefficient. Open 
squares: experimental data with the absolute volume frac- 
tion scaled to fit the Tokuyama-Oppenheim prediction [4l| . 
Open circles: same data scaled to the Beenakker-Mazur pre- 
diction [i^l. Crosses: same data scaled to maximize the over- 
lap between the experimental and numerical results for the 
structural relaxation time, r^, in the glassy regime ip > 0.55. 
The scaling factor for absolute between squares and circles 
is 1.07. 



sured ISF is w 80% or more, and we conclude that our 
DLS experiments probe essentially the self-part of the 
ISF, namely: 



/(iw,tw + t) 



(5) 



D. Determination of the volume fraction 

As mentioned above, the absolute value of the volume 
fraction after preparation of a new sample is not known 
precisely, whereas the uncertainty on the relative values 
of ip for a series of samples obtained by dilution can be re- 
duced to 10~^ or less. In most works on nearly monodis- 
perse colloidal hard spheres, the absolute scale of the 
volume fraction is obtained by preparing samples in the 
crystal-fluid coexistence region and by setting ip such that 
the experimental melting and freezing fractions coincide 
with those predicted theoretically [1, |ll| . It is important 
to remark that, unless the sample is perfectly monodis- 
perse, this method suffers from some uncertainty, since 
numerical work shows that the volume fraction at freez- 
ing depends on polydispersity. For example, in Ref. [ll| 
the authors use nominal values of Lp that are obtained by 
comparison to the freezing point of a monodisperse sus- 
pension, but they warn that the true values of tp could be 
up to a factor of 1.04 higher than the nominal ones, due 
to a sample polydispersity cr ~ 6 %. More generally, a 
relative error of about 5 % is unavoidable in experiments 
on hard sphere colloids. 
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This calibration procedure can not be applied to our 
more polydisperse suspensions, since they do not crys- 
tallize. Instead, we calibrate the absolute ip by com- 
paring the volume fraction dependence of the short-time 
self-diffusion coefficient, Ds, to theoretical predictions 
by Tokuyama and Oppenheim [4ll | and Beenakker and 
Mazur |42l |. In practice, we start from an initial guess 
of the absolute volume fractions, (/^gucss, and measure Ds 
vs. (pgucss in the range (ySguoss ^ 0.2 by fitting the ini- 
tial decay of the ISF to /(g,r ^ 0) ~ exp(— q^£>sr). 
As explained in Sec. Ill C| we measure / at qR = 4 and 
under the best index matching conditions, to make sure 
that only the self part of the ISF is probed. We then fit 
-P.sf'Pp-iiRss) to the prediction for Ds{ip) of either Ref. [4l| 
or Ref. [42|, setting cp = bipgucss- There are two fitting 
parameters in this procedure: b is the scale factor for the 
absolute volume fractions, and Dq = lim^^Q Ds{(p) the 
self-diffusion coefficient in the limit of infinite dilution. 

The results of the fit are shown in Fig. [1] together with 
the theoretical curves. The data agree well with both 
theoretical predictions, although the fit to the Beenakker- 
Mazur expression is slightly better. The ratio of the scal- 
ing factors b obtained by fitting the data to Ref. [4lj or 
Ref. [ill, respectively, is 1.07. Thus, the spread in the 
estimate of the absolute volume fractions is of order 7%, 
comparable to that in Re f. [ill . Interestingly, the present 
method was used in Ref. [43[ for a sample where calibra- 
tion using the experimental freezing point was simulta- 
neously possible, yielding fully compatible results. 

As a final, alternative procedure to calibrate the abso- 
lute volume fraction, we compare our data for the equilib- 
rium relaxation times described in more detail in Sec. IIVI 
below, to the results of Monte Carlo simulations for a 
three-dimensional binary mixture of hard spheres [2^ . 
We find that it is possible to make data in the glassy 
regime coincide over about 5 decades of relaxation times 
by adjusting the experimental scale for 1^9 by a factor very 
close to that required to fit Dg to the Beenakker-Mazur 
theory (see Fig.|T]). In this paper, as in Ref. Q, we adopt 
for convenience the scaling factor required to collapse the 
experimental and numerical relaxation times in the glassy 
regime (crosses in Fig. [1]). 

This section shows that two distinct indicators ob- 
tained by theoretical and numerical work can be used to 
adjust absolute volume fractions in moderately polydis- 
perse colloidal hard spheres, yielding results that arc con- 
sistent within a confidence interval of about Aip/ip w 7 %, 
an uncertainty comparable to that plaguing adjustments 
onto the static phase diagram due to polydispersity ef- 
fects. This uncertainty should be kept in mind when 
absolute numbers for critical volume fractions are com- 
pared between different works using different particles, 
different solvents, different techniques, and different cal- 
ibrations for the volume fraction. 
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FIG. 2: (Color online) Degree of correlation ci{t 

sured for a sample with a volume fraction ip = 0.5970. From 

top to bottom, the time delay r is 0, 1, 10, 20, 50, 100, 200, 

500, Ik, 2k, 5k, 10k, 20k, 50k, 100k, and 200k s. Note the 

logarithmic time axis, to better appreciate the initial aging 

regime. 

E. Data analysis 

As an example of the TRC data obtained in our ex- 
periments and their processing, we show in Figs [5] and 
[3] the degree of correlation c/ and the final relaxation 
of the age-dependent ISF, respectively, for a sample at 
ip = 0.5970. This is the densest system for which we 
were able to collect data at thermal equilibrium. 

The degree of correlation initially grows with t^, im- 
plying that the change in sample configuration over a 
fixed time lag r becomes progressively smaller, i.e. that 
the dynamics slows down. After about 10^ — 10^ s the 
evolution of the dynamics essentially stops and a nearly 
stationary state is reached. Figure [3] shows the two-time 
ISF obtained by averaging the c/ data of Fig. [2] over time 
windows of duration 200 s (for the younger ages) to 1000 s 
(for the oldest samples), for various ages (for clarity, 
the fast dynamics measured with the goniometer appa- 
ratus is not shown). The lines are stretched exponential 
fits to the final decay of the ISF: 

/(iw,iw + r) = A(tw)cxp{-[r/T„(tw)f (6) 

where A is the height of the plateau preceding the final, 
or a, relaxation, Tq the a relaxation time, p the stretch- 
ing exponent and i? > a small, residual base line most 
likely due to an imperfect correction of the effect of un- 
even illumination, as discussed in Ref. [s^. The aging 
behavior observed for c/ in Fig. [2] is reflected by the in- 
crease of Ta with tv,- We find that for > 5 x 10^ s the 
relaxation time does not evolve any more (see also Fig. [S] 
and the related discussion below) , while the height of the 
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FIG. 3: (Color online) Symbols: two-time intermediate scat- 
tering functions (ISFs) /(tw,tw-(-T) obtained from the degree 
of correlation shown in Fig. [2] = 0.5970). From left to 
right, tw =lk, 2k, 5k, 10k, 50k, 100k, 200k, and 300k s. The 
solid lines are stretched exponential fits to the data according 
to Eq. ®. 
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FIG. 4: Age dependence of the a relaxation time for a sample 
a.t if = 0.5970. The curves are labeled by the laser power used 
for the measurement. The decrease of Tc at large tw observed 
at the higher laser powers is an artefact due to sample heating 
and convection, which is removed by reducing the power and 
introducing a shutter. 



plateau still slightly increases, albeit extremely slowly. 
We take t^, = 5x 10"^ as the end of the aging regime and 
obtain the equilibrium value of as the average of the 
fitted relaxation time in the stationary regime. 

III. MEASURING LONG RELAXATION TIMES 

We faced several challenges due to the extremely long 
relaxation times of the system, the use of a CCD detec- 
tor, and the influence of gravity. In this section, we de- 
scribe some of the potential artefacts and the difficulties 
associated with measurements on nearly glassy samples, 
together with our solutions to overcome these problems. 

A. Sample heating and convection. 

The CCD detector is much less sensitive than a tra- 
ditional phototube or an avalanche photodiode. Con- 
sequently, a larger laser power is required, potentially 
leading to sample heating, if the particles or the solvent 
absorb light at the wavelength of the source. Under nor- 
mal gravity conditions, local heating results in convcctivc 
motion, which can significantly alter the spontaneous dy- 
namics of our samples [44| . Figure |3] shows the age de- 
pendence of the a relaxation time obtained by fitting the 
ISF to Eq. ®, for a sample at </? = 0.5970. The different 
curves are labeled by the power of the incident beam; 
when the largest available power is used (150 mW, full 
circles), Tq initially increases, as observed for many out- 
of-cquilibrium systems, but then drops sharply with in- 



creasing tvj. This surprising "rejuvenation" effect is much 
less pronounced when a reduced laser power is used (50 
mW, open squares): the slowing down of the dynam- 
ics persists over a longer time span and the decrease of 
Tq that eventually sets in is more modest. These ob- 
servations suggest that the acceleration of the dynamics, 
reported in one of our previous papers [1^, is a non- 
linear effect perturbing the equilibrium dynamics, due to 
the onset of convcctivc motion driven by sample heating. 
Indeed, we found that the solvent had became slightly 
yellowish since samples were prepared, which was respon- 
sible for increased absorption in the green region of the 
visible spectrum, where our laser operates. 

To avoid any artefact due to convection, we have added 
a shutter to the setup, so that the sample is illuminated 
only for 100 ms at each CCD acquisition, rather than 
continuously. Since the typical acquisition rate is 1 Hz or 
lower, this reduces the average power injected in the sam- 
ple by a factor of ten or more. Additionally, we perform 
our experiments at the minimum laser power compatible 
with a good beam stability (about 10 mW). Under these 
conditions, is found to increase monotonically until a 
plateau is reached, with no rejuvenation effect even after 
2 days, as shown by the triangles in Fig. [4l By changing 
the opening time of the shutter, we have checked that the 
dynamics shows no dependency on the average illuminat- 
ing power, as long as the laser power is < 10 mW and the 
shutter time is < 1 s. As a final remark on this spurious 
rejuvenation, we note that this phenomenon sets in very 
slowly, especially when the laser power is not too high. 
This hints at a subtle interplay between aging and an ex- 
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tcrnal drive, similarly to the behavior of other soft glassy 
materials to which a (modest) mechanical perturbation 
is applied during their aging [43 ■ More experiments will 
be needed to explore these analogies. 



B. Sedimentation 

Sedimentation is a potential additional source of a spu- 
rious acceleration of the dynamics, since it can act as a 
non-linear driving force. This is especially true for sys- 
tems that relax very slowly and are very sensitive to ex- 
ternal forcing [i^, \^ . In the imaging geometry of our 
CCD setup, sedimentation results in an overall drift of 
the speckle pattern on the detector, which contributes 
to the change of Ip at any given pixel and thus to the 
decay of g2. We have implemented a correction scheme, 
which will be described in detail in a forthcoming publi- 
cation; here we explain only the principles of the method. 
We use Image Correlation Vclocimetry (ICV) [i^ , a tech- 
nique similar to Particle Imaging Velocimetry, to measure 
the drift of the speckle pattern due to sedimentation for 
each pair of images taken at time and + t- The 
drift is obtained by calculating the spatial crosscorrela- 
tion between the two speckle images. If the relative dis- 
placement of the particles over the lag r is smaller than 
l/q, the lengthscale probed in a DLS experiment, the 
speckle pattern is essentially frozen, except for its overall 
drift. The second speckle image is then a shifted version 
of the first one; consequently, the crosscorrelation func- 
tion exhibits a marked peak whose position corresponds 
to the average particle drift. Once the drift is obtained, 
the second image is back-shifted numerically [i^ by the 
same amount but in the opposite direction, so as to com- 
pensate the physical drift. The corrected c/ is finally 
calculated between the first image and the back-shifted 
version of the second one. The minimum shift that can 
be reliably measured is a few hundredths of a pixel, cor- 
responding to about 50-100 nm in the sample. This is 
slightly less than 2Tr/q = 250 nm, the length scale over 
which the relative motion is probed by DLS: thus, we are 
able to detect and correct for an overall drift compara- 
ble to the relative displacement of the particles. Note 
that our correction method fails when the speckle pat- 
tern "boils" and changes faster than it drifts, since no 
peak in the crosscorrelation function can be detected re- 
liably. However, this situations corresponds precisely to 
the case where the internal, spontaneous dynamics dom- 
inates over sedimentation, making the correction for the 
drift irrelevant. 

Figure [5] shows both the raw (solid squares) and the 
corrected (open squares) 52 — 1 for the most concentrated 
sample that we have studied {(p — 0.5990) at t^, = 10^ 
s. Sedimentation clearly affects the decay of 172, which is 
about three times faster for the raw data than for the cor- 
rect ones. The inset shows the large-age behavior of the 
relaxation time extracted from raw (open squares) and 
corrected (solid squares) correlation functions. Sedimcn- 



0.8 



0.6 



+ 



^ 0.4 



^ 0.2 



0.0 



Tmr 



1 1 mill I I I mill I I I mm 1 1 1 11 



10" 



10' 



Ernrn 1 — 1 1 1 1 111 r 



: n 
n 



I L_ 



10 10 

?.(s) 







I I I 



10^ 



10 10" 



10^ 



10= 



FIG. 5: Main panel, filled squares: final relaxation of the 
ISF for a concentrated sample at <p — 0.5990 and ~ 10* 
s. The open symbols are the same data after correcting for 
the contribution of the average sedimentation velocity to the 
decay of /. Inset: for the same (fi, t„ dependence of the a 
relaxation time for the raw ISFs and the corrected data. 



tation leads to an apparent arrest of aging for tw > 10^ s, 
while the corrected data indicate that the sample barely 
equilibrates for > 10^, an age ten times larger. Note 
that, although we do correct for the average sedimenta- 
tion velocity, we can not correct for velocity fluctuations 
stemming from hydrodynamics interactions and possibly 
changing the relative position of the particles, thereby 
contributing to the loss of correlation of the scattered 
light. Since velocity fluctuations are known to be rel- 
evant in concentrated suspensions (soj . it is likely that 
sedimentation still accelerates to some extent the decay 
of (72- These data are therefore not included in the anal- 
ysis of equilibrium dynamics below. 

It is interesting to compare the potential impact of sed- 
imentation on our results to that in other works. A re- 
cent set of experiments Q on concentrated PMMA hard 
spheres has shown that aging is faster for off-buoyancy- 
matching particles than for nearly-buoyancy-matched 
colloids. Although our particles are not buoyancy- 
matched (Ap = 259 kgm~'^), they are smaller than those 
used in the optical microscopy investigation of Ref. 
enhancing the role of Brownian diffusion with respect to 
the gravity-driven drift. The relevant parameter is the 
inverse Peclet number 



A: 



47rc/Api?4 ' 



(7) 



defined as the ratio of the gravitational length to the 
particle radius (T is the absolute temperature and fee 
the Boltzmann's constant). In this work, A « 870, much 
larger than A = 44.1 for the "normal gravity" experi- 
ments of Ref. Q (for which Ap = 300 kgm~'^) and still 



8 



greater than A = 678 for their "reduced gravity" mea- 
surements (Ap = 20 kgm~^). Thus, our experiments 
are comparable to typical "reduced gravity" microscopy 
investigations. This illustrates the great advantage of 
using smaller particles to mitigate gravity effects, due to 
the R^^ dependence of A. In order to further investi- 
gate the role of gravity in the slow dynamics of colloidal 
hard spheres in the large A regime, light scattering mea- 
surements on small particles that are partially buoyancy- 
matched or microgravity experiments will be needed. 

As a final remark on potential artefacts, one may won- 
der whether mechanical instabilities could be partially 
responsible for the ultra-slow decay observed in Fig. [5l 
where the relaxation time can be as large as 1.5 days. 
To test this possibility, we have measured g2 for a sam- 
ple compressed at random close packing and where very 
little dynamics is expected to occur. Although 172 does 
exhibit some decay, whose origin is still unclear (either 
mechanical instability or some extremely slow rearrange- 
ment of the packed particles, or the contributions of a 
small number of "rattlers"), we stress that this relax- 
ation is much slower than that observed for the most 
concentrated sample studied in this work. Furthermore, 
the relaxation of the most concentrated sample (Fig. [5]) is 
a factor of 10 slower than the slowest equilibrium relax- 
ation reported here and in Ref . Q : we can therefore rule 
out any significant artefact due to mechanical instability. 
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FIG. 6; Age dependence of the a relaxation time for samples 
at volume fractions 0.5876 < f < 0.5990, as indicated by 
the labels. Different symbol fillings correspond to indepen- 
dent experiments on the same sample. Large-fw data for the 
two largest volume fractions are corrected for sedimentation 
eff'ects. The line has a slope of one, showing that the initial 
growth of Ta is in general close to linear. 



C. Non-equilibrium aging effects 

Besides artefacts due to sample heating and sedimenta- 
tion, the dynamics of very concentrated samples become 
intrinsically difficult to measure, because of aging and dy- 
namical heterogeneity. Aging is a quite general feature 
of glassy systems (see e.g. [5l| and references therein), 
where the dynamics slows down with time, or sample 
"age", <wj as the system evolves towards equilibrium. In 
traditional DLS, the intensity correlation function has 
to be extensively averaged over time, making it impos- 
sible to capture the evolution of the dynamics for non- 
stationary systems. By contrast, the TRC approach and 
other multispeckle methods [13, HI, \EM allow one to fully 
characterize time-evolving dynamics, because the time 
average is replaced in part, or completely, by an average 
over the slightly different q vectors associated to distinct 
pixels of the CCD. 

Figure shows the age dependence of for samples 
prepared at different volume fractions (p > 0.5876. Sym- 
bols with the same shape but different filling correspond 
to independent experiments at the same volume fraction. 
Data for the two largest ip are corrected for sedimentation 
effects at the largest t^,; in all other cases, corrected and 
raw data yield the same results. Initially, increases 
with tv,; in this regime, the growth of vs. tw is gener- 
ally close to linear, as shown by the line oc iw For the 
most diluted sample in Fig. [S] (92 = 0.5876), the initial 
aging regime is very short and is comparable to the time 



needed by the sample to equilibrate at the temperature 
imposed by the copper holder; consequently, no precise 
aging law can be determined. After the initial regime, 
Ta saturates and becomes almost independent of time, 
suggesting that all samples reach equilibrium, with the 
possible exception of the most concentrated suspension 
{(fi = 0.5990, squares in Fig. for which the dynam- 
ics seem to keep slowing down over the full duration of 
the experiment (more than 3 days), although at an in- 
creasingly smaller rate. Note that the experiments are 
quite generally well reproducible, with the exception of 
the earliest stages, where presumably the dynamics is 
more influenced by the exact initial configuration and 
thermalization. 

The time needed for reaching equilibrium, teq, is about 
30 times greater than the (equilibrium) relaxation time 
for the samples at the lowest volume fraction, (p = 0.5876. 
However, this values is certainly overestimated due to the 
contribution of the thermalization time, of the order of 
10^ s. For the samples at ip ^ 0.5953 and 0.5970, for 
which full equilibration is reached and where the initial 
thermalization time is negligible compared to all relevant 
time scales, teq/Ta ~ 3 — 4. This confirms the intuitive 
notion that equilibrium is reached on the time scale of 
a few structural relaxation times, i.e. after a few cage 
escape processes. Note that we are able to show that 
there is indeed an equilibrium regime even for samples 
whose a relaxation time is more than 7 orders of mag- 
nitude larger than that of very dilute suspensions. This 
significantly extends the range of samples that do reach 
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FIG. 7: Age dependence of the stretching exponent p obtained 
by fitting the final decay of the ISF to Eq. ([SJ . Same symbols 
as in Fig. [6] 

equilibrium, compared to previous work [ll[ , where sam- 
ples with Ta larger than 10^ times the Brownian decay 
time were classified as aging samples, most likely due 
to technical difficulties to access much larger relaxation 
timcscalcs. 

The time evolution of the stretching exponent p{t^,) is 
shown in fig. [71 For sufficiently large , p saturates at 
0.5-0.6, regardless of the volume fraction of the particles 
and the history of the samples. This stretched expo- 
nential shape might be related to sample polydispcrsity, 
since in previous experimental work on similar but less 
polydispersc colloidal hard spheres (cr w 6%) the ISF 
was shown to relax nearly exponentially [T]] |. However, 
we point out that we find a similar value in the simula- 
tions described below in Sec. IIVBI for both a = 5.77% 
{p fa 0.55) and a = 11.5% {p ~ 0.65). Interestingly, for 
the most concentrated sample p remains well above the 
asymptotic value measured for all other samples, further 
suggesting that this sample has not fully relaxed during 
the experimental time window. Note that here we never 
observe "compressed" exponential relaxations (p > 1), in 
contrast to the data of Ref. [s^ . We recall that for the 
latter laser heating yielded convection, which was most 
likely responsible for faster-than-exponcntial relaxations. 



D. Dynamic heterogeneity 

A detailed discussion of dynamical heterogeneity in 
samples at intermediate volume fractions (0.20 < ip < 
0.576) has been presented in Refs. [^,[5^, [5^. Here, we 
simply report that both the degree of correlation c/ and 
the agc-depcndcnt ISF become more erratic and exhibit 
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FIG. 8: (Color online) Degree of correlation for the most con- 
centrated sample {ip = 0.5990). Note the large fluctuations of 
ci, especially for tv, < 5 x 10"^ sec, indicative of a temporally 
heterogeneous dynamics. From top to bottom, r is 0, 5, 10, 
25, 50, 100, 250, Ik, 2.5k, 5k, 10k, 25k, 50k, 100k, and 150k 
sec. 

significant temporal fluctuations as ip grows. An example 
is shown in Fig. [8] for c/ , and in Fig. [9] for the correspon- 
dent timc-resolvcd ISFs. The data presented here are 
for the most concentrated sample that we have studied, 
(f = 0.5990. Note that the fluctuations of cj are particu- 
larly relevant at early ages, when the sample is presum- 
ably farther from its equilibrium configuration. Because 
of these fluctuations, the shape of f{t^„,t^ + T) is less well 
defined, and the fits are generally poorer. However, we 
stress that Tq remains reasonably well defined, since the 
decay of the ISF, although somehow "bumpy", is gener- 
ally not too stretched, with the exception of the earliest 
ages. 

The observation that c/ and / arc increasingly noisy 
as increases is consistent with the recently reported 
growth of cooperatively rearranging regions in super- 
cooled colloidal suspensions approaching the glass tran- 
sition m, [13, [53I, in analogy with the behavior of 
molecular glass formers [H, . Indeed, as the size of the 
regions that rearrange cooperatively increases, the num- 
ber of such regions in the scattering volume decreases: 
the dynamics are averaged over a smaller number of sta- 
tistically independent objects, leading to enhanced fiuc- 
tuations [6l|. The fact that dynamical fluctuations at 
high ip are significant suggests that cooperatively rear- 
ranging regions may extend over a sizeable fraction of 
the scattering volume. This would imply a correlation 
length of the dynamics of the order of several hundreds 
of particle diameters, much larger than the ~ 10 parti- 
cles reported in previous works [H, 113, Ull at lower 
(p, but comparable to recent findings in polydisperse col- 
loids near random close packing |62| . More experiments 
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FIG. 9: (Color online) Two-time ISFs obtained from the 
cj data shown in Fig. [8] (v? = 0.5990). From left to right, 
i„ = 500, 5k, 7k, 10k, 20k, 30k, 50k, 70k, 100k, and 170k s. 
Dynamical heterogeneity at early times results in noisy ISFs, 
whose shape significantly departs from stretched exponential. 
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FIG. 10: (Color online) Inset: representative ISFs for equi- 
librated samples at volume fractions 0.0480, 0.3096, 0.4967, 
0.5555, 0.5772, 0.5818, 0.5852, 0.5916, 0.5957, 0.5970, from 
left to right. The lines are stretched exponential fits to the fi- 
nal decay of /. Main figure: scaling of the a relaxation of the 
ISFs for ip > 0.4967 when using the reduced variable t/tci. 
The line is the function 0.68exp[-(r/T,,)''-5^**]. 



will be required to confirm these intriguing preliminary 
results. 



IV. ANALYSIS OF EQUILIBRIUM DYNAMICS 
A. Equilibrium relaxation times 

The first two sections were meant to convince the 
reader that it is possible to access the equilibrium dynam- 
ics of our colloidal hard sphere sample over a broad range 
of timescales for several well-controlled volume fractions. 
We now turn to the analysis of these results, which were 
briefly presented in a shorter paper Q. 

In the inset of Fig. [10] we show the time decay of the 
ISF at selected volume fractions from a dilute system at 
If ~ 0.05, up to very large volume fractions, ip = 0.597 
(same data as in Ref @). Our data cover a broad range of 
11 decades in time, and we follow the slowing down of the 
equilibrium dynamics over about 7 decades in relaxation 
times. 

In agreement with previous work [l^ . we find that 
time correlation functions decay exponentially when vol- 
ume fraction is moderate, with a time constant that in- 
creases weakly with (p. When (p is increased above some 
'onset' volume fraction, ip w 0.5, the relaxation becomes 
strongly non-exponential, with a two-step decay that be- 
come more pronounced as p continues to grow, signalling 
the arrest of particle motion for intermediate timescales. 
While the short-time decay is relatively unaffected by the 
increase of p>^ the characteristic time of the second decay 



corresponding to the structural relaxation of the fluid 
increases rapidly over a small range of volume fraction. 
For (p > 0.597, we were not able to collect data for which 
stationary behavior and thermal equilibrium could be un- 
ambiguously established, and this volume fraction marks 
thus the location of the 'experimental colloidal glass tran- 
sition', pg « 0.6, by analogy with the temperature scale 
Tg in thermal glasses. We do not attach any particular 
signiflcance to pg, because it obviously depends on the 
details of the experimental set-up and of the sample. 

As described in Sec. Ill El we fit the long-time decay 
of the time correlation functions to a stretched exponen- 
tial form, Eq. The resulting fits are shown in the 
inset of Fig. [TO] as continuous lines. They describe the 
data very well and yield quantitative confirmation of the 
above qualitative remarks. We find that the amplitude 
A of the decay is A w 1 at low ip when relaxation is 
mono-exponential, so that p w 1 in that regime. When 
p increases, we find that A drops to A 0.7, signalling 
a two-step process. At the same time, the stretching ex- 
ponent p decreases and stabilizes around p ~ 0.56, indi- 
cating that structural relaxation occurs through a broad 
distribution of timescales. 

The fact that p is only weakly dependent on the volume 
fraction suggests that the structural decay of the ISFs 
should collapse when plotted against rescaled time, in the 
spirit of the 'time temperature superposition principle' 
frequently observed in molecular glass-formers. We show 
such a scaling in the main plot of Fig. [TOl where the ISFs 
for p > 0.4967 are plotted as a function of r/ra, with 
the structural relaxation time issued from the stretched 
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FIG. 11: (Color online) Open squares, left axis: volume frac- 
tion dependence of the a relaxation timescale for equilibrated 
samples. The black line is an MCT fit to the data, Eq. ([l}. 
Open circles, right axis: a relaxation timescale measured by 
van Megen et al. in [ll[, in units of the Brownian time, 
Ti, = i?r /(6Do). The blue line is an MCT fit to the data 
of [l3|. Crosses: same data as the circles, plotted as a func- 
tion of volume fraction scaled by a factor of 1.031, a value 
smaller than the uncertainty on if> stated in [ll|]. The scaled 
data satisfactorily overlap with the measurement of this work. 



exponential fit to the data. The data collapse reasonably 
well onto a single master curve, confirming that the time 
volume fraction superposition principle holds for colloidal 
hard spheres, as observed in [ll|. 

We present in Fig. [11] the evolution of the equilibrium 
relaxation time Tq,(</3) extracted from the fits of the ISFs. 
For the sake of clarity, only data for > 0.48 are shown: 
data at all volume fractions were reported in Ref. Q. 
While Ta increases by about 1 decade from ip = 0.05 
to (fi ~ 0.5 (data not shown), the structural relaxation 
slows down dramatically by about 6 decades when ap- 
proaches the largest volume fraction considered in this 
work. Following previous analysis [Tl1 |. we fit Tq, to an 
algebraic divergence, as in Eq. ([T]), using ipc and 7 as 
adjustable parameters. Such a power law, predicted by 
MCT, cannot account for our data over the whole range 
of volume fraction, and the range of volume fractions fit- 
ted by Eq. ^ must be chosen with great care. As is also 
found in molecular glass-formers, we find that the power 
law Eq. ([1]) describes a window of relaxation times of 
about three decades, immediately following the onset of 
glassy dynamics, as shown in Fig. [TlJ The values of the 
exponent, 7 = 2.5 ± 0.1, and of the critical volume frac- 
tion, ific ~ 0.590 ± 0.005, depend on the precise range of 
volume fraction fitted, but they agree well with previous 
work [ll[ , and theoretical and numerical analysis 0, [IJ] , 
as we will further discuss in Sec. IIVCI If we attempt to 
include the largest volume fractions in our fit, we find 
that 7 takes unphysically large values, up to 7 w 6. We 
interpret this finding as the sign that the growth of Tq, is 
not appropriately described by Eq. ([1]) at large volume 



fraction. 

A more direct indication can be observed in Figs. [TO] 
and 1111 since we were able to collect equilibrium data at 
volume fractions ip above the fitted value of the critical 
volume fraction cpc- Therefore, the dynamic singularity 
implied by Eq. ([1]) is not observed in our sample, and 
ergodic behavior can be detected for ip > ipc- Just as is 
gencrically found in thermal glasses, we conclude there- 
fore that the singularity predicted by mode-coupling the- 
ory is avoided in our experimental colloidal system. Al- 
though suggested by theory and computer simulations, 
such an observation was not reported in experimental 
work before. 

As discussed in detail in Ref. , we find that an expo- 
nential divergence of Tq, accounts for our data very well, 
as in Eq. . We argued in Ref. Q that the best fit of our 
data was obtained for C = 9.8 x 10~^, Tqo = 6.5 x 10^^ 
sec, 6 = 2.0± 0.2 and po = 0.637 ± 0.002, although 
imposing 6 = 1, as in free volume predictions, yields 
an acceptable fit and a smaller critical volume fraction, 
Pq{S = 1) = 0.614 ± 0.002. 

Interestingly, the non-trivial exponent 5 w 2 is also 
supported by the numerical data for the polydisperse 
system of quasi-hard spheres studied in Sec. IIVBI (see 
Fig- [H] below), and by Monte Carlo results for a binary 
mixture of hard spheres [§, [l^l (also shown in Fig. [T2)) . 
Additionally, a recent scaling analysis of the glass tran- 
sition occurring in systems of soft repulsive particles 
yields S = 2.2 ± 0.2 for the hard sphere limit [2J, [6^, a 
value compatible with that obtained here from the fit of 
Eq. Thus, there is mounting evidence from both sim- 
ulations and experiments that a simple VFT law, Eq. ^ 
with (5 = 1, is not the most accurate description of the 
dynamics of hard spheres. 

Overall, we find that while the onset of dynamical slow- 
ing can be described by an MCT divergence at a critical 
volume fraction pc, upon further compression a crossover 
from an algebraic to an exponential divergence at a much 
larger volume fraction pQ is observed, showing that the 
apparent singularity at pc does not correspond to a gen- 
uine 'colloidal glass transition'. Instead the system en- 
ters a regime where dynamics is 'activated', in the sense 
that Ta increases exponentially fast with the distance to 
Pq. These observations therefore suggest that the glass 
transition scenario in colloidal hard spheres strongly re- 
sembles the one observed in molecular glass-forming ma- 
terials. 



B. Polydispersity effects 

A central finding from our experiments is the absence 
of a genuine dynamic singularity occuring at the fitted 
mode-coupling singularity pc- Since polydispersity influ- 
ences the dynamics of colloidal hard spheres, we must ask 
how our findings may change for samples with a different 
polydispersity. 

In Ref. Ill a sample less polydisperse than ours (<t w 
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6% as opposed to cr « 10%) was analyzed using simi- 
lar light scattering techniques, but no deviation from an 
algebraic singularity was reported. Additionally, exper- 
iments performed on samples at different polydispersi- 
ties showed that increasing the polydispersity at constant 
volume fraction produces an acceleration of the dynam- 
ics [131 • These observations suggest a possible explana- 
tion for our findings. One could argue that the observed 
ergodic behavior in our sample at large (f stems from the 
fact that the glass state above the critical volume frac- 
tion determined in Ref. pjj is 'melted' by polydisper- 
sity effects. In other words, polydispersity could push 
ific to a much larger value so that our fluid states could 
all be observed below the "true" characterizing our 
sample. This view is apparently supported by numerical 
analysis [35j , which showed that volume fractions for the 
location of random close packing could be shifted from 
Vrcp = 0.64 to (^icp = 0.67, from monodisperse to 10 % 
polydisperse samples, respectively, thus suggesting a sim- 
ilarly large shift in the MCT critical volume fraction. 

This string of arguments is immediately contradicted 
by the fact that we have self-consistently determined ipc 
for our own sample, and that we do observe ergodic be- 
havior for if which are larger than this critical volume 
fraction. It remains thus to establish that polydispersity 
does not qualitatively affect the dynamics in such a dras- 
tic manner that the results reported in this article do not 
shed light on the behavior of less polydisperse samples. 

To definitely settle this issue, we performed computer 
simulations of a quasi-hard sphere system at different 
polydispersities [6^. Following previous work [6^, we 
study an assembly of point particles interacting via a 
purely repulsive potential: 



2R^ 



36 



(8) 



where r^j is the distance between particles i and j, and 
Rij ~ {Ri + Rj)/2 with Ri the radius of particle i. The 
prcfactor e is an energy scale. The interaction poten- 
tial is very steep, and should therefore constitute a good 
approximation of the hard sphere potential. We set the 
polydispersity by drawing the particle radii from a flat 
distribution: Ri e [R_- SR/2, R + SR/2], so that the 
average diameter is Ri = R and the polydispersity is 



cr = - R2/R = 5R/{VV2R). We have studied the 

system for two values of polydispersity, 5R/R = 0.2 
and 0.4, leading to polydispersities a ~ 5.77 % and 
a ~ 11.5 %, which compare well with the colloidal sam- 
ples studied in Ref. [ll| and in the present study, respec- 
tively. 

We follow previous work dedicated to the Hamilto- 
nian ([5]) in the context of hard spheres. We solve 
Newton's equations for N particles of mass m inside a 
three-dimensional simulation box of size L with periodic 
boundary conditions [6^ . For the inverse power law po- 
tential in Eq. ([5]), density and temperature can be com- 
bined in a unique control parameter, F = (/sT"^/^^, so we 
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FIG. 12: (Color online) Top panel: self-intermediate scat- 
tering function from numerical simulations of a quasi-hard 
sphere system at two different values of the polydispersity, a, 
and various volume fractions. For the less polydisperse sys- 
tem, if = 0.053, 0.50. 0.55, 0.58 and 0.592, from left to right. 
For the more polydisperse system, Lp — 0.074, 0.50, 0.55, 0.58, 
0.5914 and 0.60, from left to right. Bottom panel : relaxation 
time Ta for the two systems, with data from Monte Carlo sim- 
ulations of a binary mixture superimposed. The three sets 
of data are fitted to the MCT prediction, Eq.(IT]), using the 
same exponent 7 = 2.7. The critical volume fraction issued 
from the fit is shown in the label. 



fix the temperature and energy scales. ksT = e = 1/3, 
and vary the system size to change the volume fraction, 
if = 47riVi?3[l-H(5i?/4i?)2]/(3L3). Since we deal with soft 
spheres, it is of course not possible to compare absolute 
values for critical volume fractions with results obtained 
in the true hard sphere limit. Timescales are expressed 
in units of 2Ry^m/e. We choose the self- intermediate 
scattering function Fs{q, t) as a dynamic observable, and 
work at a single wavevector, qR = 3.9, close to the first 
diffraction peak. 

We perform the simulations in two steps. We pre- 
pare initial configurations at the desired state points, and 
perform a long equilibration run using periodic velocity 
rescaling to reach thermal equilibrium. We then perform 
a production run in the microcanonical ensemble. We 
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check very carefully that thermal equilibrium is indeed 
reached during equilibration. We also check that our re- 
sults do not depend on the duration of these two steps, 
and we repeat them for at least five independent samples 
at each volume fraction to obtain better statistics. The 
main issue, mentioned in the original work of Ref. [66j , is 
the possibility for the less polydispersc system (cr w 6%) 
to crystallize in the course of the simulation. In this case 
the sample is discarded and a new sample studied. We 
found that this problem is much more severe than re- 
ported in Ref. [6y], and so we were not able to collect 
data above = 0.587. For the more polydisperse sys- 
tem (ct « 12%), we never detected crystallization and 
we were able to collect data up to 1^9 = 0.60, where we 
stopped because we could not reach equilibrium within 
our numerical time window. It is interesting to note that 
crystallization similarly prevented the colloidal samples 
of Ref. [ll| to be studied at very large volume fractions, 
while no such problem is encountered for the colloidal 
system studied in this article. 

We report in the top panel of Fig. [T^] some representa- 
tive results for the self-intermediate scattering function 
of the system at two polydispersities, and different vol- 
ume fractions. This figure is qualitatively very similar 
to the experimental data of Fig. IIOI with a pronounced 
non-exponential relaxation which slows down dramati- 
cally when ip increases. When volume fraction is low, we 
see that the change in polydispersity has no detectable 
influence on the dynamics, and data for the two systems 
nicely superimpose. A first visible effect of polydisper- 
sity is that we do not have data for the 6 % polydisperse 
sample at very large values of the relaxation time because 
the system eventually crystallizes. A second effect can be 
detected when the glassy regime is entered. For a fixed 
(p, the dynamics of the system become faster when poly- 
dispersity is increased, as observed previously [63. m^. 

To quantify these observations, we fit the long time 
decay of the relaxation to a stretched exponential form, 
as in Eq. ([6|) , and show the results for Tq {(f) in the bottom 
panel of Fig.[T21 As mentioned before, the polydispersity 
has very little effect when dynamics is fast, Lp < 0.5, while 
the more polydispersc sample relaxes faster in the glassy 
regime. 

We have fitted the dynamic slowing down with increas- 
ing Lp to the algebraic form suggested by MCT. As ex- 
pected, we find that polydispersity affects the value of 
the critical volume fraction, although the change is rather 
modest, since ipc increases from (pc ~ 0.595 to (pc ~ 0.597 
when changing polydispersity from 6 % to 12 %. Note 
that we used the same exponent 7 = 2.7 to fit both sets 
of data over a comparable window of relaxation times, 
showing that polydispersity does not quantitatively af- 
fect the form of the slowing down. It must be noted that 
the shift in the absolute value of (pc is much less severe 
than expected from that for the random close packing 
volume fraction on the basis of the numerical work of 
Ref. m. 

It is interesting to note that deviations from the MCT 



power law cannot be detected very clearly for the less 
polydisperse system, because crystallization prevents ob- 
servation of larger relaxation timescales. For the more 
polydisperse system, however, we find that deviations 
from the MCT power law at even larger volume fractions 
become observable, since a broader window of relaxation 
timescales can be measured. We believe that similar ef- 
fects are relevant also in experiments. As we will show 
shortly, both the experiments reported in Ref. [ll[ and 
our data are described by a comparable power law over a 
similar window of relaxation timescales. However, while 
previous work did not detect deviations from an algebraic 
divergence, by studying a more polydisperse system, we 
can efficiently suppress crystallization and access a range 
of equilibrium relaxation timescales that are beyond the 
volume fraction regime which can be described by MCT, 
but the MCT regime itself is not qualitatively affected 
by polydispersity effects. 

For completeness, we report in the bottom panel of 
Fig. [T^] the Monte Carlo results obtained in Ref. Q for a 
three-dimensional 50:50 binary mixture of hard spheres of 
diameter ratio 1.4, for which polydispersity is cr « 16.7 %. 
Again, we find that the data is qualitatively unaffected 
by the use of a very different particle size distribution, by 
the fact that a true hard sphere potential is used, and by 
the use of a stochastic microscopic dynamics. An MCT 
power law also describes the data over a comparable time 
window of about three decades, with a similar exponent 
7 = 2.7, and a critical volume fraction ipc = 0.593. Re- 
member that we cannot directly compare the result for 
ipc to the ones for the quasi-hard sphere system, since 
the latter is dependent on the temperature scale used 
(a larger temperature would yield a smaller critical den- 
sity). Overall, it is reassuring that whenever comparison 
is possible, neither the choice of a microscopic dynam- 
ics, nor the particle size distribution and polydispersity 
seem to influence the dynamics in a drastic manner. This 
certainly suggests that quantitative comparison between 
different experimental samples is meaningful. 



C. Comparison with previous work 

Considering colloidal hard spheres from the viewpoint 
of the glass transition field, it is very natural to expect, 
just as is found for all molecular glasses, that the dy- 
namic singularity deduced from fitting data to an alge- 
braic law predicted by MCT is eventually avoided as the 
glassy regime is entered more deeply. Yet, the present ex- 
periment is the first direct demonstration that the MCT 
dynamic transition is avoided in a colloidal hard sphere 
sample. Although we are tempted to assume that our 
conclusions should apply to all colloidal hard spheres 
samples, this statement would seem to require additional 
discussion. 

First, it is reassuring to observe that several com- 
putational studies performed in recent years with hard 
spheres, using different types of polydispersities, and dif- 
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ferent types of microscopic dynamics indicated the pres- 
ence of strong deviations from an algebraic divergence 
when approaching Lpc- In that sense, the experiments of 
Rcf. Ill| stood as an exception, although arguably an 
important one. 

In Fig.[TTl we report data from Ref. [ll| {t'^, blue cir- 
cles and crosses and right axis), together with our data, 
taken from Ref. (ra, black squares and left axis). The 
open circles are the data as presented in Ref. [ll|, using 
<f values that may be underestimated by a factor < 1.04, 
because of polydispcrsity, as discussed in Sec. HID I An 
MCT fit to Tl^i^p) yields (^c = 0.572 = ipf^ ^^{hhie 
line). The crosses are the same data plotted against vol- 
ume fraction scaled by a factor ^pc/^p^ = 1.031, with 
Lpc = 0.59 the critical volume fraction obtained by fitting 
our data to the MCT prediction, Eq. ([T]) [t^ . Note that 
the scaling factor required for the MCT critical volume 
fractions to coincide is smaller than the uncertainty in the 
absolute Lp for both sets of data. When compared over 
the range of volume fraction where MCT applies, both 
sets of data are nearly indistinguishable, confirming that 
in this regime sample polydispcrsity has a very small in- 
fluence on the dynamics, as indicated by the simulations 
discussed in Sec. IIVBI 

Figure [TT] shows also that the range of tinicscalcs re- 
ported in Ref. [llj is relatively smaller by about two 
decades than the present data. Assuming that the be- 
havior of the fluid at larger Lp for the sample of Ref. [ll| 
would also compare well with the one we observe, one 
realizes that the data point of Ref. [ll| at the highest tp 
lies just below the volume fraction where deviations from 
MCT predictions should start becoming visible, possibly 
explaining why no such deviations were observed in that 
work. Of course, obtaining equilibrium data at larger 
Lp for this sample would be difficult, because slow dy- 
namics would start competing with crystallization, as we 
observed numerically in Sec. IIVBI for the sample at the 
lowest polydispcrsity. 



V. CONCLUSIONS 

We have discussed some of the challenges that must 
be faced when investigating the slow dynamics of con- 
centrated colloidal systems, together with the solutions 
we devised. These challenges include the determination 
of the absolute volume fraction and avoiding or at least 
minimizing artefacts due to sedimentation and convec- 
tion. At very large volume fractions, aging and dynami- 
cal heterogeneity make measurements even more difficult. 
The data become more erratic, suggesting that a proper 
ensemble average is not being achieved. 

The data presented here and in |9i] show unambigu- 
ously that, for our sample, the power-law divergence of 
Ta{ip) is avoided and that at high (p the growth of the 
relaxation time is well described by a generalized VFT 
law. A comparison with the measurements reported by 
van Megen et al. [ll| indicates that, in the MCT regime. 



the two sets of data are compatible within the experi- 
mental uncertainty on the absolute volume fraction. This 
suggests that the slow dynamics of hard spheres is only 
weakly sensitive to the precise value of polydispcrsity, at 
least up to moderate values a ^ 10%. Our numerical 
work supports this conclusion. The agreement between 
our data and previous work strongly suggests that the 
absence of a true algebraic divergence of Tq should be a 
general feature in colloidal hard spheres, independent of 
the details of the system. In this respect, hard spheres 
should be similar to molecular glass formers, for which 
the MCT transition is avoided as well. 

Several questions remain open. While we were able 
to rule out the existence of an MCT glass transition at 
p = ipc ^ 0.59, our experiments do not allow us to 
determine unambiguously whether the location of the 
divergence of Tq, obtained by a generalized VFT fit, 
(/So = 0.637, is distinct from (/Srcp- Addressing this is- 
sue would allow one to discriminate between competing 
theoretical approaches [6§|, e.g. recent excluded volume 
approaches [l6l and thermodynamic glass transition the- 
ories [itI Il8l|. However, providing an experimental an- 
swer will be particulary difficult: on one hand, ipo can 
only be obtained by extrapolation. On the other hand, 
(ySrcp is difficult to locate both on a conceptual stand- 
point [2^ and because any deviations from the behavior 
of an ideal hard sphere potential may become important 
very close to (p-cc-p- 

Another promising line of research concerns the aging 
regime. We have shown that at the highest volume frac- 
tions investigated the dynamics slows down over a period 
of the order of a few (equilibrium) structural relaxation 
times. Whether a similar behavior also applies to more 
diluted samples remains to be ascertained. Indeed, in our 
experiments the thermalization time, during which con- 
vective motion is likely to set in due to small temperature 
gradients across the sample, was too long to reach any 
conclusive result for (p < 0.5953. The large dynamical 
heterogeneity of the most concentrated sample suggests 
that dynamical fluctuations may be much larger than 
previously thought in out-of-equilibrium samples at very 
high ip. The interplay between dynamical heterogeneity 
and the evolution of the dynamics in the aging regime 
appears as a promising field for future investigation. 
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